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Motivated by recent experimental data on thin film superconductors and oxide interfaces we pro- 
pose a random-resistor network apt to describe the occurrence of a metal-superconductor transition 
in a two-dimensional electron system with disorder on the mesoscopic scale. We explore the in- 
terplay between the statistical distribution of local critical temperatures and the occurrence of a 
lower-dimensional (e.g., fractal-like) structure of a superconducting cluster embedded in the two- 
dimensional network. The thermal evolution of the resistivity is determined by an exact calculation 
and, for comparison, a mean- field approach called effective medium theory (EMT). Our calculations 
reveal the relevance of the distribution of critical temperatures for clusters with low connectivity. 
In addition, we show that the presence of spatial correlations requires a modification of standard 
EMT to give qualitative agreement with the exact results. 
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I. INTRODUCTION 

Although two-dimensional (2D) electron systems have 
always been a topical subject in condensed matter 
physics, further renewed interest in this field has been 
triggered by the recent discovery of superconductivity 
at the metallic interface between two insulating oxide 
layersp™^ In addition, experiments on certain thin con- 
ventional superconducting filmP^ reveal pseudogap ef- 
fects, a phenomenon which could unveil new features of 
the superconducting transition in low-dimensional disor- 
dered systemsP In the particularly interesting cases of 
homogeneously disordered superconducting titanium ni- 
tride or indium oxyde thin films, scanning tunneling spec- 
troscopy data display inhomogeneities in the local den- 
sity of states on mesoscopic spatial scales, clearly indicat- 
ing the existence of an inhomogeneous superconducting 
state. Since tunneling experiments cannot directly access 
the 2D metallic layers at the oxide interfaces (although 
some attempts made on samples with very few top layers 
have detected interesting inhomogeneous textures^) mag- 
netization and transport measurements are the only way 
to investigate these systems. Also in this case, evidences 
for an electronic phase separation have been found in 
LaA10 3 /SrTi0 3 layers Moreover, the sheet resistance 
curves ifcl(T) near T c , both in LaA10 3 /SrTi0 3 [2] and 
LaTi0 3 /SrTi0 3 [3] systems, show a peculiar feature of a 
marked tail on the low-temperature side. 

Despite their differences, oxide interfaces and thin films 
share the common feature of a broad superconducting 
transition. In a previous workEQwe investigated the pos- 
sible origin of this pronounced width and found that su- 
perconducting fluctuations (alone) cannot account for its 
occurrence, while a model of quenched mesoscopic disor- 
der can well produce the observed broad transitions. In 
this previous work we also found that "tailish" features 
can occur in the resistivity curves when disorder has spa- 
tial correlations. Specifically, while uncorrelated disorder 



can generically account for the pronounced width of the 
metal-superconductor transition, tails only occur when 
disorder displays a correlated character or when it oc- 
curs on nearly one-dimensional subsets of the system. 
It is therefore important to understand how (correlated) 
disorder can act as to give rise to the tailish resistivity be- 
cause this feature can be a signature of specific physical 
mechanisms at work in these systems. One natural pos- 
sibility is that disorder correlations and/or nearly one- 
dimensional (or fractal) structures arise in the system 
from nanofractures, topological defects like crystal dis- 
locations, steps, twinning domains, and so on. Another 
possibility is that low-dimensional superconducting sub- 
sets can spontaneously arise in the electronic gas even 
in the presence of a homogeneously distributed disorder. 
In particular, one should remember that, besides fluctu- 
ations of the energy gap A, disorder can cause the local- 
ization of electrons, transforming an otherwise metallic 
system into an insulator. If the metal is also a super- 
conductor then, at low temperatures, dis order can in- 
duce a superconductor-insulator transition! 12 ! 13 ! In this 
case, a theoretical proposal has been put forward, which 
predicts that superconductivity at high disorder is main- 
tained by a fragile coherence within a small set of pre- 
formed Cooper pairs that are characterized by an anoma- 
lously large binding energy. Consequently, in the vicinity 
of the superconductor-insulator transition, both the in- 
sulator and the superconductor contain these preformed 
Cooper pairs that either localize, leading to an insulating 
state, or condense into a coherent state with a glassy be- 
havior around the transit ion E2HH1 The superconducting 
state then occurs on a small (possibly fractal) subset of 
the system, which forms in the electron gas around T c . 

To shed light on the physics of disordered low- 
dimensional superconductors, we carry out a systematic 
investigation of a 2D metallic system with quenched dis- 
order represented by a random- resistor network (RRN). 
As in our previous work of Ref. , the resistors repre- 
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sent mesoscopic metallic islands where superconductivity 
occurs at a local critical temperature T c . The supercon- 
ducting regions are assumed to be large enough to have a 
fully established local coherence and to make the charg- 
ing effects negligible. As a consequence, a mutual phase 
coherence immediately establishes between two neighbor- 
ing superconducting islands, as soon as both have become 
superconducting. This clearly distinguishes our frame- 
work from the case of granular superconductors, where 
the grains have usually nanoscopic sizes of the order of 
the coherence length of the pure system. In this way we 
disregard phase-fluctuation effects in order to unambigu- 
ously identify the effects of disorder (and of its specific 
character). In particular, to understand the effects of 
spatial correlations, we investigate here the effects of two 
fractal-like distributions of the disordered resistors. In 
this way, while completely disregarding the underlying 
mechanisms leading to such correlated disordered struc- 
tures, we focus on their phenomenology. 

The paper is structured as follows. In Section II we in- 
troduce the generic aspects of our model, while in Section 
III we report the results of our numerical calculations for 
various kinds of fractal structures and of disorder distri- 
butions. Section IV contains a discussion of the results 
and our concluding remarks. 

II. THE RANDOM-RESISTOR MODEL AND 
EFFECTIVE MEDIUM THEORY 

In a previous worlP^ we considered a model of a 2D 
electron gas with mesoscopic defects as a square lattice 
whose bonds are assigned a random resistivity pi. More 
precisely, each bond is assigned a local superconducting 

(i) 

transition temperature T c , extracted from a given dis- 
tribution W(T C ), and the resistivity of the bond is written 
as pi — po@(T — Tc : ), with the same high-temperature 
value po on all bonds, 0(x) being the Heaviside step func- 
tion. By decreasing the temperature T, more and more 
bonds become superconducting, and global superconduc- 
tivity establishes as soon as a percolating superconduct- 
ing path is formed in the system. 

The simplest descript ion of a RRN is provided by the 
effective medium theorjEEO (EMT), which is a mean- 
field treatment replacing the random resistors pi with an 
effective medium resistivity p em such that 

y Pi -Pern =Q> (1) 
t-r 1 Otpi + p em 

where the parameter a is related to the connectivity of 
the network. For a cubic network in D spatial dimen- 
sions, a = D — 1. 

In the previous work of Ref. [TT] (to which we refer the 
reader for further detail), we used EMT as a benchmark 
and compared p em with the exact numerical determina- 
tion of the resistance of the RRN, obtained solving the 
Kirchhoff equations for the network in the presence of a 



difference of potential V between two opposite sides of 
a N x N square lattice. Once the current I(T) flowing 
through the network at a temperature T is obtained, the 
resistance is found as R(T) = V/I(T). We showed that 
EMT provides a very good description of the tempera- 
ture dependence of the resistivity, independently of the 
distribution of T c , provided disorder is spatially uncorre- 
cted. 

T he sol ution of Eq. (ffl), specialized to our case, 
reads™! " 

p em {T) = (l + a)p 9(T-T a ) jf dT c W{T c ). (2) 

Here, T a is the critical temperature of the effective 
medium and is determined by the equation 

w s (T a )= dT c W{T c ) = — — , (3) 

Jr a 1 + a 

where w s {T) = /+ °° dT c W(T C ) is the statistical weight 
of the superconducting bonds at a temperature T, mea- 
suring the frequency of occurrence of bonds with T c > T. 
From Eq. (pi) it is readily seen that p em —>■ po for T — > oo 
and p em {Ta) = 0. 

In the present paper we consider a model in which the 
superconducting bonds only form a spatially correlated 
subset of the whole system. This subset provides a "skele- 
ton" network of bonds which become superconducting 
below a random local T c , while the embedding system 
is a metal with temperature-independent resistive bonds. 
We shall explore the joint effects of the connectivity of 
the network and of the statistical distribution of the local 
critical temperature W(T C ). We shall mainly rely on the 
exact numerical solution of the Kirchhoff equations, but 
shall refer to the EMT results to provide an interpreta- 
tion to the outcomes of our calculations. When dealing 
with EMT, the presence of spatial correlations requires 
the introduction of an effective connectivity a, different 
from the standard connectivity. 

III. EFFECTS OF FRACTAL-LIKE 
DISORDERED STRUCTURES 

A. The general framework 

To systematically investigate the effects of correlated 
disorder, we consider the occurrence of a cluster where 
the zero-resistance state is carried by only a minor set 
of superconducting regions, reducing the effective dimen- 
sion of the 2D electron gas. We disregard the physical ori- 
gin (classical, as due to mechanical stresses, non-uniform 
growth and so on, or quantum-mechanical, as - for ex- 
ample - associated with the occurrence of coherence on 
a thin cluster of Cooper pairs near the superconductor- 
insulator transition) of this low-dimensional supercon- 
ducting "skeleton" and we translate it into the framework 
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of the RRN. Specifically, we consider a RRN on which 
fractal- like structures are superimposed, i.e., clusters ex- 
hibiting strong spatial correlations and scale invariancc 
(although the scale invariance is limited, once the struc- 
ture is implemented on a finite - e.g., 100 x 100 - lattice). 
Each bond of these spatially correlated clusters is given 
a critical temperature extracted from a probability dis- 
tribution W(T C ), while the bonds not belonging to the 
"fractal" cluster have fixed resistances po (throughout the 
paper we take po = 1). By varying the geometry of the 
clusters, as well as the probability distribution, the influ- 
ence of these two aspects on the temperature dependence 
of the resistivity is studied. Note that in this approach 
the spatial distribution of superconducting bonds is con- 
sidered to be independent from the probability distribu- 
tion of the critical temperatures. This is a simplifying 
assumption allowing for a more general and systematic 
analysis of the separate effects of geometry and of sta- 
tistical distribution. This assumption should of course 
be dropped whenever the physical mechanisms selecting 
the "fractal" superconducting cluster also determine its 
T c distribution. 

It is worth pointing out that in the previously consid- 
ered models of uniform uncorrelated disorder 11 , global 
superconductivity occurred as a true percolative phe- 
nomenon: Upon lowering the temperature, more and 
more bonds became superconducting until a percolative 
cluster of superconducting bonds was reached. The ge- 
ometrical support and the general features of this tran- 
sition were typical of the percolation transition. On the 
other hand, in the present model the geometric support 
of bonds which can become superconducting is provided 
by a predetermined fractal-like subset of the network. 
When T is lowered, more and more bonds on this fractal 
support become superconducting and global supercon- 
ductivity only occurs when a percolating path through 
the whole system is formed inside the subset. In this 
case, for instance, the fractal dimension of the supercon- 
ducting support is not due to pure percolation, but it 
clearly depends on the dimensionality of the supporting 
fractal-like subset. This "percolation inside a fractal" is 
what gives rise to global superconductivity in the present 
model. 

In order to cover a wide range of possible T c distribu- 
tions we choose the two extreme cases of Gaussian and 
Cauchy statistics. The choice is motivated by their differ- 
ent asymptotic behavior: While the first has small tails 
and well defined momenta at all orders, the second is 
characterized by strong tails causing all momenta except 
the mean to diverge. We take a Gaussian defined by the 
mean value p\ and variance a, 

W(T C ) = _^ e -{T.-^?/2*> (4) 
V27T<7 

and a Cauchy (i.e., Lorentzian) distribution defined by 
the mean value pi2 and the width 7, 

Wm " .[V + ( r.- w )-| (5) 



Throughout the paper the same parameters p\ = pi = 1, 
a = 0.1 and 7 = 0.08 will be adopted. The parameters a 
and 7 are chosen so that the resistivity curves obtained 
with the Gaussian and Cauchy distribution have the same 
slopes close to the transition when space correlations are 
absent (Fig. [I). 
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FIG. 1. (Color online) Resistance curves of the uniform clus- 
ter. The blue curves correspond to a Gaussian distribution 
with fi = l,cr = 0.1 and the green to a Cauchy distribution 
with y, = 1,7 = 0.08. The markers report the thermal evo- 
lution of the resistance obtained through the exact numerical 
calculations and the full lines the EMT solution. 

As recalled in Sec. |nj one of the findings of Ref. [TTj 
was that in the absence of space correlations EMT per- 
forms remarkably well in reproducing the exact results of 
the RRN, irrespective of the specific T c distribution. As 
an illustration, we report in Fig. [I] the resistivity curves 
of a uniform system obtained with the distributions here 
considered: The Gaussian (blue curves) and the Cauchy 
distribution (green curves). One notices that the resis- 
tance curves obtained with EMT match very well the 
exact numerical solutions. Thus, in the absence of cor- 
relations, EMT provides a useful tool to investigate the 
p(T) dependence. In particular, the theory states that 
the resistivity of a RRN vanishes as soon as the weight 
of superconducting bonds w s exceeds 1/(1 + a), where a 
is the connectivity of the system. In the present case of 
a uniform two dimensional lattice, a = 1 and the perco- 
lation threshold at the transition yields w s = 1/2. Re- 
ferring to the resistance curves in Fig. [T] one notices 
that the critical temperature of the system is T a = 1, 
coinciding with the mean value of the (symmetric) dis- 
tributions where half of the bonds have become super- 
conducting. Thus, we correctly reproduce the predicted 
value of w s = 1/2. At high temperature, the heavier tail 
of the Cauchy distribution results in a greater number 
of bonds with switched-off resistance. The resistivity for 
the Cauchy distribution remains smaller until the tem- 
perature T a = 1 is reached where w 1 ^ — wf = 1/2, the 
superscripts C, G referring henceforth to the Cauchy and 
Gaussian distribution, respectively. 

To give a quantitative picture of the resistivity close to 
T a one can take the derivative of Eq. ^ defined in the 



4 




framework of EMT: 

p' em (T) = (l + a)p W(T c = T), (6) 

for T > T a , which shows that the slope of the resistivity 
is proportional to the value of the distribution at equal 
temperatures. Combining this equation with the previ- 
ous result w s (T a ) — 1/2, where the distribution has a 
maximum, one concludes that it is impossible to obtain 
"tails" in the present case. In fact, due to the bell-shape 
of the distributions the slope of the p(T) is maximal at 
the mean value and thus also at the critical tempera- 
ture of the system. In the same line of reasoning, one 
realizes that a tailish behavior can be obtained even for 
a spatially uncorrelated distribution, in the very specific 
case of a symmetric bimodal distribution of T c 's (see Ref. 
[TT]). Here, we disregard this possibility and focus in- 
stead on the issue of space correlations. 

B. Diffusion Limited Aggregation 

The first fractal-like cluster we implement is obtained 
through a simple growth process, generated by Brown- 
ian motion and known as diffusion limited aggregation 
(DLA). Its construction is quite simple: A particle is re- 
leased at the left edge of a 2D lattice and let diffuse to 
the right. More precisely, the particle moves one bond 
to the right and then with equal probability one bond 
up or down. This procedure is iterated until the parti- 
cle stops, as soon as it reaches the top, bottom or right 
edge, where it sticks. Then, other particles are launched 
in sequence and halted either when reaching one of the 
three edges or a bond already occupied by one of the 
previously diffused particles. The cluster obtained in a 
250 x 250 square lattice after diffusing 50000 particles is 
defined by the bonds where the particles sticked. Due to 
a saturation at the left edge, the total number of super- 
conducting bonds only amounts to about 25000. Once 
this large cluster is obtained, we select a 100 x 100 sub- 
lattice as shown in Fig. [2] and perform the calculations 
for this smaller cluster. In this way, we try to model 
a more physical case where the fractal covers the whole 
sample. So henceforth, we consider a 100 x 100 lattice, 
where only bonds belonging to the cluster are assigned 
a critical temperature T c . The other bonds form a resis- 
tive background and are assigned the resistivity po at all 
temperatures. 

Fig. [3] reports the resistance curves calculated on the 
100 x 100 DLA cluster. We remark that for T > 1 
the curves for the Gaussian and Cauchy distribution are 
rather similar. Lowering the temperature, the system 
starts to exhibit a stronger dependence on the specific 
statistics. In the Gaussian case, a percolating path is al- 
ready formed around T a = 0.78, whereas for the Cauchy 
distribution one has to go as low as T a = 0.33 for the sys- 
tem to become fully superconducting. Examining the for- 
mation of the superconducting cluster, one finds that at 
T a , the fraction of superconducting bonds are wf = 0.99 



FIG. 2. Spatial distribution of superconducting bonds ob- 
tained through diffusion limited aggregation. The red square 
indicates the 100 x 100 lattice which is used for the calcula- 
tions. 



and w]? — 0.97 respectively. So practically all the bonds 
have to be superconducting in order for the phase tran- 
sition to occur. 

The reason for this strong condition lies in the very 
small set of percolating paths, i.e., the effective connec- 
tivity of the cluster is very close to zero. In other words, 
the cluster has a rather marked one-dimensional charac- 
ter and few missing superconducting bonds are enough 
to prevent the system from percolating. To give a quanti- 
tative estimate of the effective connectivity of the cluster 
we define the following measure (inspired by Eq. ([3| de- 
fined in the EMT framework), 



For instance, comparing (T a ) — 0.5 and 
wf' DLA (T a ) = 0.99, where the superscripts U and 
DLA refer to the case of a uniform and DLA super- 
conducting cluster, respectively, one realizes that the 
statistical weight of the superconducting bonds clearly 
depends on the geometry of the cluster and, in particular 
close to the critical temperature T a , on its effective 
connectivity In addition, we shall show below that 
taking a e J[ A = 1/0.99 - 1 = 0.01 the EMT restricted 
to the cluster produces p(T) curves close to the exact 
results. Therefore, even though space correlations do 
not enter explicitly, equation Q still provides a good 
estimate for an effective connectivity, once the relevant 
space correlations are accounted for by restricting EMT 
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to the cluster. A second effect of the small set of 
percolation paths is a strong dependence on the specific 
realization of the distribution. Since the number of 
bonds in the percolating cluster is at most of order N D 
(where D is the fractal dimension of the cluster), two 
different realizations of the same distribution might be 
quite different. This effect is particularly severe for the 
Cauchy distribution, which has larger deviations from 
the mean. Since the Gaussian distribution, having small 
deviations from the mean, produces a superconducting 
cluster with an effective connectivity a e J[ A — 0.01, we 
assume this value as representative of the appropriate 
connectivity of the underlying DLA cluster. On the 
other hand, the realization of the Cauchy distribution 
of critical temperatures on the cluster in Fig. [2] is such 
that a percolating path forms at finite temperature 
but, calculations of p(T) for different realizations of 
W(T C ) show that on average the phase transition does 
not occur. This result can be understood by noticing 
that the Cauchy distribution has rather wide tails so 
that wf(T < 0) 0.025. Thus the fraction of bonds, 
which never become superconducting is substantial 
and, owing to the low connectivity of the DLA cluster 
(requiring wf LA — 0.99), it may well happen that 
non-superconducting bonds prevent full percolation to 
be established down to zero temperature. 
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FIG. 3. (Color online) Resistance curves of diffusion limited 
aggregation cluster. The blue curves correspond to a Gaussian 
distribution with p = 1, a = 0.1 and the green to a Cauchy 
distribution with p = 1,7 = 0.08. The markers report the 
thermal evolution of the resistance obtained by exact numer- 
ical calculations, the dashed lines by EMT on the system and 
the full lines by EMT on the cluster. 

This effect makes it quite apparent that on correlated 
clusters, the low values of connectivity (like ct e J[ A = 
0.01) renders the system very sensitive to the statistics 
of the critical temperatures. 

Fig. [3] also shows that the EMT, when applied to 
the whole RRN (which is formed both by the fractal- 
like cluster and its complementary background), clearly 
misses the correct temperature dependence of the resis- 
tivity. According to its construction, outlined in Sec. 



p"H this mean-field-like theory neglects all spatial correla- 
tions and this discrepancy is expected. In addition, we 
remark that the decrease of p with T is a finite-size effect: 
As soon as we introduce spatial correlations, the cluster 
dimension D becomes smaller than 2, and the statisti- 
cal weight of superconducting bonds is smaller (or equal) 
than N D /N 2 which tends to zero for A — > 00. 

One can try to circumvent the shortcomings of EMT by 
evaluating it on the bonds belonging to the DLA cluster 
only. With this adjustment, the relevant space correla- 
tions are effectively incorporated and EMT with the ef- 
fective connectivity reproduces quite accurately the main 
features of the exact solution: The regular behavior for 
T > p and the stronger dependence on the distribution at 
low temperatures, due to the cluster's low connectivity. 
However, there is an overall shift from the restricted EMT 
solution with respect to the exact calculation. At high 
temperatures, where only very few bonds have become 
superconducting, it is not favorable to force the current 
to flow through many resistive bonds simply to reach a 
few superconducting bonds. Therefore, the current still 
flows through the system in parallel (only locally per- 
turbed by the few superconducting bonds present). By 
decreasing the temperature the superconducting cluster 
increases and consequently, also the amount of current 
going through it. Close to the transition the shift then 
nearly vanishes because the resistive background ceases 
to play an important role. More precisely, the current is 
essentially carried by one or a few large quasi-percolating 
paths, which only need few resistors to switch off in or- 
der to fully percolate. For these paths, it is immaterial 
whether they are embedded in a two-dimensional (mark- 
ers in Fig. [3]) or simply a quasi one-dimensional system 
(full lines in Fig. [3). 

C. Symmetrized Random Walk 

The topology of the DLA cluster considered in the pre- 
vious section is characterized by few backbones (forming 
connected paths between the two vertical edges of the 
lattice) and many dangling branches. Now, we explore a 
different situation in which the fractal-like geometry does 
not contain any dangling branches. 

We modify the construction scheme of the DLA in the 
following way: Instead of keeping only the final position 
of the diffused particles, their entire trajectory is incor- 
porated into the cluster. In addition, the particles are 
free to cross the paths of previously diffused particles. 
Launching 50 particles from each of the four edges and 
letting them diffuse perpendicular to the initial edge, a 
cluster of the form shown in Fig. [4] is obtained. 

As it was the case for the DLA, the resistance curves of 
the symmetrized random walk (SRW) shown in Fig. [5]are 
very similar for T > 0.9, the small difference originating 
simply from the difference in the distributions. Look- 
ing at the superconducting cluster formed at T = 0.9, 
one finds that the weights of superconducting bonds are 
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FIG. 4. Spatial distribution of superconducting bonds ob- 
tained through a symmetrized random walk on a 100 x 100 
lattice. 




FIG. 5. (Color online) Main panel: Resistance curves of sym- 
metrized random walk cluster. The blue curves correspond 
to a Gaussian distribution with p — 1, a = 0.1 and the green 
to a Cauchy distribution with fi = 1,7 = 0.08. The markers 
report the thermal evolution of the resistance obtained by ex- 
act numerical calculations, the dashed lines by EMT on the 
system and the full lines by EMT on the cluster. The inset 
shows a comparison of the exact resistivity curves obtained 
with a Cauchy distribution (p — 1,7 = 0.08) for the uniform 
(circles), the DLA (pluses) and the SRW (squares) cluster, 
respectively. 
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= 0.84 and 



,,c 



0.70. respectively. Interestingly, 



this difference has no effect on the resistivity. In fact, at 
this temperature, the size of connected superconducting 
regions is of the order of 10 bonds for both distributions. 
In this regime, having a few superconducting bonds more 
or less does not produce a noticeable difference in the re- 
sistivity. 

This behavior clearly changes as one approaches the 
transition temperature. For T slightly larger than T a the 
current is carried by long-range superconducting regions 
which lack only few bonds to fully percolate. Inspection 
of the superconducting cluster growth reveals that the 
percolation threshold is reached as soon as wf(T^) — 
0.92 and wf{T%) = 0.90. The closeness of these two 
values is again an indication that the threshold weight 
w s (T a ) only depends on the "geometry" (i.e. its effective 
connectivity) of the underlying cluster and not on the 
specific T c distribution. While the threshold is reached at 
— 0.86 in the first case, the heavy tail of the Cauchy 
distribution obliges the system to go as low as =0.76 
in the second case. 

These considerations on , w^ LA and wf RW sug- 
gest that the high-temperature behavior is mainly de- 
termined by the dimensionality of the cluster. In fact, 
coarse-graining analysis using the box-counting method-^ 
shown in Fig. ^reveals that the dimension of the clusters 
is D DLA sa jjShw ^ 18 This ig renec ted in the congru- 
ence of p DLA {T) and p SRW (T) for T > 0.9, as well as in 
their leftward shift with respect to p u (T) (with D u = 2), 
as is shown in the inset of Fig. [5] 

In contrast, the calculation of the connectivity at the 
transition point reveals that OL e J^ w « 0.1, which corre- 
sponds to a nearly one-dimensional cluster. So, despite 



its 2D appearance, the low-temperature region, being 
more sensitive to the effective connectivity than to the 
dimensionality, behaves as quasi-one-dimensional. As far 
as the comparison between the EMT and the numerical 
exact results is concerned, the conclusions are analogous 
to the discussion given for the DLA cluster. 

To develop a quantitative understanding of p(T) close 
to the transition, we turn to an analytic approach based 
on a coarse-grained picture of the cluster. We consider 
the random walk cluster as a uniform system made up 
by large conducting segments containing M ^> 1 bonds. 
In the framework of EMT the percolation threshold is 
reached when 1/2 of the segments are superconducting. 
One can reformulate the condition in terms of the single 
bonds; in order for a large segment to become supercon- 
ducting, all its single bonds have to be superconducting, 



i.e. 



P M = 1/2 



(8) 



where p is the probability for a single bond to be su- 
perconducting. Close to the global critical point, the 
resistivity can be approximated by 



1 



V 



(9) 



with C > 0. The purpose of the model is to analyze the 
possible occurrence of a tail in the resistivity curve, i.e., 
the behavior of the slope of p(T) at T w T c : 



dp 
dT 



TaT r 



-Mp M - l p' 



(10) 
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FIG. 6. (Color online) Determination of fractal dimension 
using the box-counting method: The lattice is subdivided into 
N box = (N/L box ) 2 boxes of size A box = L box x L box . For each 
box k, we calculate the ratio \ih = N bond /N c , between the 
number of bonds belonging to the portion of the cluster inside 
this box N bond and the total number of bonds belonging to the 
cluster N c . The quantity R q = Ylk=T Mfc ^ s then related to 



the box size according to R q 



rD(q-l) 
J b o X 



where D is the fractal 



dimension. The figure shows a log-log plot with L box varying 
from 1 to N/2 (circles) . The full lines are linear regressions of 
In R q= 2 oc D In L box ; the DLA cluster (green curve) and the 
SRW cluster (blue curve) have both slope D = 1.8. 



where Cpo was set to 1 for simplicity and p 1 — dp/dT. 
Since T w T c and M 3> 1 the following approximations 
can be made, 



1 
2 

exp 
M = 



1 

ln(2) 



ln(2) 



1 ln(2) 

M y ' 



1 



V 



(11) 
(12) 

(13) 



Inserting the two expressions into Eq. |10[ we get an ex- 
pression for the derivative of the resistivity as a function 
of the single bond probability close to the critical point, 



dp 
dT 



T«T C 



2 K, p(l-p) 



■iln(2)-iL_ 
2 'l-p 



-ln(2)- 
2 V ; 1 



W(T C 



J£ W(T c )dT c 



(14) 
(15) 
(16) 



In terms of the distribution the condition p ss 1 (see 



Eq. (14)) means that the leading behavior is obtained 
taking the limit W(T C — > — oo) (or, more precisely, p. — 
T c a, 7; we point out that the critical temperature 
stays finite) . This implies that the slope of the resistivity 
is proportional to the ratio of two small numbers. For the 



Gaussian distribution, the ratio is 



W(T C ) 



/2-n-er 



1-J-W(f c) df c l-(l-^) 
M \T C - Mil 



y/2a 



2a 2 



(17) 



where the short-hand notation x — (T c — p,\)/\J 7 2.o was 
used. For the Cauchy distribution, the ratio is 



W(T C ) 



^r[7 2 + (Te-M2) 2 ] 



l-/~W(f c )df c l-(l-^) 



1 



\Tc - H2\ 



(18) 



One observes that the slope for the Gaussian distribution 
depends linearly on the deviation from the critical tem- 
perature, while in the case of the Cauchy distribution the 
dependence is inversely proportional. The evaluation of 
the expressions p7| and ( [l8| with the results of the SRW 
yields p'(T c G ) sa p[ T c) ~ 0.15, which is smaller than the 
exact results given in Fig. |5j The reason for this mis- 
match lies in the assumption M 3> 1, while in the present 
case M w 5 . Nevertheless, the calculations highlight the 
strong influence of the distribution for low-dimensional 
structures with M 1. 

To investigate the SRW and illustrate the analytic ap- 
proach presented above, we briefly consider a cluster ob- 
tained through coarse graining. The cluster contains ev- 
ery fifth column and fifth row of the original network. 
In this way, one obtains a uniform grid equivalent to the 
uniform case of a iV/5 x N/5 lattice where each bond is 
made up by 5 bonds of the original network. 




FIG. 7. (Color online) Resistance curves of coarse grained 
cluster. The blue curves correspond to a Gaussian distribu- 
tion with jj, = 1, a — 0.1 and the green to a Cauchy distribu- 
tion with fi — 1,7 = 0.08. The markers report the thermal 
evolution of the resistance obtained by exact numerical calcu- 
lations, the dashed lines by EMT on the system and the full 
lines by EMT restricted to the cluster. 



Comparison of Fig. [7] and Fig. [5] reveals a marked re- 
semblance between the p{T) curves for the random walk 
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and the coarse grained cluster. This shows that the reg- 
ular coarse-grained cluster indeed captures the main fea- 
tures of the strongly disordered SRW cluster. 



IV. CONCLUSIONS 

Our analysis is purely phenomenological in character: 
We embed a given fractal-like structure in a purely metal- 
lic environment, we impose a given T c distribution on 
it, and then we calculate the resulting resistivity curve. 
While, as we already mentioned, this approach lacks any 
(possible, but not mandatory) microscopic connection 
between the geometrical and the T c distribution, we sep- 
arately access the distinct effects of these two ingredients. 

As a first result, we would like to point out that the 
geometrically correlated character of the superconduct- 
ing clusters considered here greatly degrades the perfor- 
mance of the EMT mean-field-like approach. According 
to our experience in spatially uncorrelated disordei^Q, 
this failure has nothing to do with the strength of dis- 
order or the relative density of the superconducting and 
the non-superconducting bonds, but is a mere result of 
space correlations. On the other hand, we showed that 
it is possible to find a good qualitative agreement by re- 
stricting EMT to the cluster. In this case, one needs to 
introduce an effective connectivity (see Eq. |7])). 

Calculating the standard (local) connectivity given by 
a = {z)/2 — 1, where (z) is the average number of near- 
est neighbors, we obtain aoLA = 0.38 and asttw — 0.15. 
The comparison with a e J[ A — 0.01 and a e /^ w = 0.1 
reveals that the effective connectivity is very different 
from the locally defined connectivity. This is quite natu- 
ral since we are considering transport phenomena, which 
are ruled by paths connecting distant regions of the sys- 
tem. In other words, the presence of regions with large 
internal connectivity, is rather immaterial for transport 
if this regions are weakly connected to one another. 

We also find that the effective connectivity is an intrin- 
sic property of the cluster and is independent of the T c 
distribution. This result is quite natural because upon re- 
ducing T, the T c distribution only determines how "fast" 
the resistances are switched off. However, this is irrel- 
evant as far as the amount of superconducting bonds 
needed to percolate is concerned. 

The main outcome of our analysis is that the effects of 
space correlations render the superconducting transition 
strongly dependent on the disorder distribution, and in 
particular on its low-temperature asymptotics. This is 
not surprising: As soon as percolation happens through 
a quasi one-dimensional path, the system needs to ex- 
plore (almost) the entire T c distribution, and thus be- 
comes very sensitive to its low-temperature part. This 
effect can be made particularly clear within the analytic 



approach of Sec. Ill C where the slope of the resistivity 
at T c was shown to depend in inverse ways on the dis- 
tance of T c from the mean of the Cauchy or Gaussian 
distribution. 

After the detailed analysis of the various effects of 
space correlations and T c statistics we come back to the 
initial question of whether and which of these aspects is 
relevant in explaining the tailish behavior of p(T) in oxide 
interfaces. Based on the above results a possible ingredi- 
ent is indeed a fractal-like structure of superconducting 
bonds. Due to its high dimensionality (i.e. larger than 
a purely one-dimensional cluster), the cluster consists of 
enough superconducting bonds to cause a linear decrease 
of the resistance with a relatively high slope. Close to 
the percolation threshold, however, the crucial property 
determining the shape of the transition is not the dimen- 
sionality but the effective connectivity. There, the quasi- 
one-dimensional behavior of the DLA and SRW clusters 
results in a more or less pronounced tail of p(T > T a ). 
The shape of the tail is then controlled (to some extent) 
by the low-temperature asymptotics of the distribution. 

Our analysis and approach are therefore a useful tool 
to explore the physical mechanisms at work in the super- 
conducting oxide interfaces, where a precise fitting of the 
resistance curves provides indications on both the intra- 
grain mechanisms leading to a local T c and to the inter- 
grain connection ruling the establishing of the global su- 
perconducting state via the percolative "chaining" of su- 
perconducting regions on the fractal support (see Sec. 
Ill A for the specific meaning of the word "percolative" 
in the present model). The intragrain superconductiv- 
ity manifests itself in the high-temperature region of the 
resistivity, where the restivity starts to bend down when 
isolated superconducting regions are created. In this case 
the T c distribution and the relation between its width 
and the high temperature resistivity might provide in- 
formations on the intra-grain pair formation (Cooper 
pairs in the presence of quenched impurities,™ disor- 
dered bosoni c pref ormed pairsj^ glassy superconduct- 
ing transition! 14 ! 15 *) On the other hand, near the global 
T c , the second percolative phenomenon is ruled by the 
geometric character (in particular by the effective con- 
nectivity) of the superconducting cluster, and might be 
informative about the origin of the inhomogenous metal- 
superconductor structure with coexisting metallic and 
superconducting islands. 
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